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Abstract 

We discuss the ground state and some excited states of the half- 
filled Hubbard model defined on an open chain with L sites, where only 
one of the boundary sites has a different value of chemical potential. 
We consider the case when the boundary site has a negative chemical 
potential — p and the Hubbard coupling U is positive. By an analytic 
method we show that when p is larger than the transfer integral some 
of the ground-state solutions of the Bethe ansatz equations become 
complex- valued. It follows that there is a "surface phase transition" 
at some critical value p c ; when p < p c all the charge excitations have 
the gap for the half-filled band, while there exists a massless charge 
mode when p > p c . 
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The Mott-insulator transition is a fundamental phenomenon where the 
strong correlation among electrons plays an essential role. The existence of 
the insulating phase, which we call the Mott insulator, can not be explained 
within the standard framework of the band theory. For the ID Hubbard 
model, it is well known that under the periodic boundary condition, the 
charge gap exists only for the positive Hubbard coupling U > and at 
half-filling. [ l| Near the transition point, however, the system shows quite 
nontrivial many-body effects. @, [| f|, §] For instance, the effective mass 
diverges at half-filling for the Hubbard ring. 

In order to investigate many-body effects near the transition point very 
precisely, let us consider a Hubbard chain in which only one site has a differ- 
ent chemical potential. With the local chemical potential we can effectively 
change the number of electrons (or holes) of the Mott insulator, infinitesi- 
mally. Let us assume L electrons in the Hubbard chain with L sites. The 
system is divided into two parts; a "surface" part consisting of only the site 
with the local chemical potential, and a "bulk" part of the other L — 1 sites. 
When the local potential is zero, the number of electrons in the bulk part is 
given by L — 1; when it is very large, no electron should occupy the surface 
site and hence all the electrons should be in the bulk part. Thus, by control- 
ling the parameter, the effective number of electrons in the bulk part (L — 1 
sites) can be changed continuously from L — 1 to L. The property of the 
electrons in the bulk part is unique: if we consider a standard closed system, 
the electron number will be given by some integer and can not increase or 
decrease infinitesimally. 

In this paper, we consider the Hubbard system defined on an open chain, 
where one of the two boundary sites is chosen as the surface site. We discuss 
how the half-filled ground state changes under the local chemical potential. 
We derive complex ground-state solutions of the Bethe ansatz equations by 
an analytic approach. We find them explicitly for some finite-size systems, 
solving the Bethe ansatz equations numerically. Then, we calculate the en- 
ergy of the ground state with the complex solutions, analytically. It is our 
hope that the study of this paper might shed some light on some new aspects 
of the many-body effects of the Hubbard system near the metal-insulator 
transition. 

The study of this paper could be related to some real ID systems such as 
Cu-0 chain |J and quantum or atomic wires |7|, ||. The open-boundary ID 
Hubbard system with the boundary chemical potential could be realized in 
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some ID Hubbard system in reality, where the local chemical potential may 
play the role of a nonmagnetic impurity or a bias potential. 

Let us introduce the ID Hubbard Hamiltonian under the open-boundary 
condition, in which only the 1st site has the local chemical potential —p. 
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3=1 <r=U i =1 cr=U 

Here Cj i<7 and rij jC stand for the annihilation and number operators of electron 
located at the jth site with spin a, respectively. We recall that U denotes 
the Hubbard interaction and t the transfer integral. Hereafter we set t = 1. 
The Bethe ansatz equations for the ID Hubbard model have been discussed 
under some different cases of open boundary conditions. || fL(| [0], [T^| (See 



also Ref. JT3|) In this paper we discuss the open-boundary Hubbard system 
with p > 0. 

For iV electrons with M down spins, the roots of the Bethe ansatz equa- 
tions are given by momenta (charge rapidities) kj for j=l to N and rapidities 
(spin rapidities) v m for m — 1 to M. With some functions Z^{k) and Z s L {y\ 
the Bethe ansatz equations can be written as 

Z c L (k J )=I j /L for j = l,...N, Z s L (v m ) = J m /L for v = l,...,M. 

(2) 

Here the quantum numbers Ij and J m are given by some integers. 

Let us consider the half-filled band under the boundary chemical poten- 
tial, where N = L and M = L/2. Hereafter we assume that L is even. 
We consider analytic continuations of the functions Z1(k) and Z s L [v) with 
respect to the parameter p. Let us introduce an adiabatic hypothesis that 
the quantum numbers Ij and J m should be constant when we continuously 
change the parameter p. Under the hypothesis, all the solutions of the Bethe 
ansatz equations can be labelled by their quantum numbers. When p = 0, we 
can order the ground-state roots k/s and t> m 's such that Ij — j 'for 1 < j < L 
and J m = m for 1 < m < L/2. The hypothesis is consistent with our ana- 
lytic arguments and numerical results. Thus, for any value of p, the quantum 
numbers of momentum kj and rapidity v m are given by j and m, respectively. 

We now consider the Bethe ansatz equations more explicitly. Let us 
denote by I max (/ m » n ) the largest (smallest) integer of the quantum numbers 
of real momenta over all possible excitations and by J max (Jmin) that of real 
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rapidities. Then, the set A£ e (A£ e ) of all the possible quantum numbers I/s 
(J m 's) for real momenta (rapidities) are given by 

^re {I mini ^min 1; • • • 5 Ifnax} (^re {'Jmini Jmin 1; • • • 5 Jmax}^) ■ (3) 

Let us write by A c hole (A^ oie ) the set of the quantum numbers of holes of 
real momenta (rapidities) in the ground state. Then, the set of the quantum 
numbers of the real momenta (rapidities) for the ground state is given by 
A< = A c re - A c hole (A* = A* e - A s hole ). Let us denote by P max {J max ) the 
largest integer of the set A^ (A*). Then, I^ ax < I ma x, in general. We 
introduce the symbol A c im (Af m ) for the set of the quantum numbers for 
complex-valued momenta (rapidities) in the ground state. In terms of the 
symbols, the functions Z1{k) and Z s L {v) for the ground state are written as 
follows 

2k 1 1 

Zl(k) = ^ + Z E E 1 (smk-rv n ) + j-z c B (k), 

n€A| r=±l 

Z s L (v) = yE E 9 1 (v-rsmk j )-j E E - ru n ) + y^(u), 

^ jeAg r=±l ^ neA| r=±l ^ 

(4) 

where the functions z#(/c) and 2#(i>) are given by 

z |(jfc) = — ^ log f 1+petk ) + y y eAsmk-rVm), 

z s B ( v ) = E E 0i(u-rsinA;j) - E ^^{v-rv m ). 

(5) 

Here, the functions 9 n (x) have been defined by 9 n (x) = 2 tan 1 (x/ (nu)) / (2ir), 
where u is given by u — U/A. An outline of the derivation of the Bethe ansatz 
equations is given in Appendix A. 

When p is larger than some critical values of p, some of the ground-state 
solutions become complex-valued. The number of complex roots is different 
for four regions of p, which are divided by the critical values p c /s. They are 
given by p c \ = 1, p C 2 = u + \J\ + u 2 , p c ^ = 2u + VT+ Au 2 . Let us introduce 
some notation. We define symbol k by k — log \p\ for p > and p < 0. 
We also define a by a = sinh/t/w for p > and p < 0. The notation of 



the critical points is summarized as p C j = (j - l)u + Jl + (j - l) 2 u 2 for 
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j = 1,2,3. If a set A is empty, we denote it by A = 0. Then, the sets of 
quantum numbers are given by the following. 

1. For < p < p c i, we have no boundary solutions. The sets of quantum 
numbers are given by 

A^ = {1,2,...,L}, A* = {l,2,...,L/2}, 

A c = A s = AS , = A? , —6. 

im im hole hole V" 

The Jmin's are given by the following 

Imin 1, Imax Li J min 1, Jmax L/1. 

2. For pd < p < p c2 (0 < a < 1), we have a complex- valued momentum 
k L given by 

= 7T + IK — i§L- (6) 

The sets of quantum numbers are given by 

A° = {2,...,L}, A* = {l,2,...,L/2}, 

AL = {1}, AL = AU = Kole = <P- 
The I m inS are given by the following 

Imin 1, Imax L 1, Jrnin 1, Jmax L/1. 

3. For p c2 < p < p C 3 (1 < a < 2), we have the complex momentum ki 
and the complex rapidity V\ given by 

ki = n + iK — iSL, 

Vi = i(a — l)u + irji. (7) 

The sets of quantum numbers are given by 

A^ = {2,...,L}, A* = {2,...,L/2}, 

A? m = {L}, AL = {1}, AU = AU = 0. 
The Jmin's are given by the following. 

Imin 1) Imax 1, 'Jrnin 2, Jmax L/1. 
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4. For p C 3 < p (2 < a), we have the following three complex roots ki, ki 
and v\ 



k x = ihgl K (a-2)u + ^/(a-2) 2 u 2 + 1) + iS 1} 

k L = TT + iK-i5 L , (8) 

vi = i(a — l)u + irjx. 

We call them a boundary k — A string. The sets of quantum numbers 
are given by 

A° = {2,3,...,L-1}, A' = {2,3,...,L/2}, 

A^ = {1,L}, AL = {1}, A^ e = {L}, A^ e = 0. 
The I m inS are given by the following 

Imin 2, Imax -^- / ; J rain 2, Jmax L/2. 

We note that when p > p c3 , a hole of real momenta appears in the half-filled 
ground state at I — L\ l max = L — l and I max = L when p > p c3 . We also note 
that 5i,8l and r\\ are exponentially small except for some neighborhoods of 
the critical points. For instance, we can show b\ = 0(p~ 2L ) for p cl < p < p c2 . 
The quantities Sl, Si and rji are explicitly evaluated in Appendix B. 

For the case when p < 0, some complex boundary solutions have been 
discussed for the ID Hubbard model under the open-boundary conditions 
[|T|, [15], |I6| , where the quantum numbers of the complex rapidities k 2 and 
v x correspond to 1\ = 1, J2 = 2 and J\ = 1, respectively. Furthermore, when 
the band- width t is very large and the electron density A^/L is very small, 
the boundary solutions fc l5 A; 2 and f x for the case of p < can correspond to 
the boundary solutions of the ID interacting spin-1/2 Fermi system, which 
had been discussed in Ref. |17|]. (See also Appendix B.) 

For the half-filling case, the ground-state energy for p > is related 
to that of p < through the particle-hole transformation, which will be 
discussed in Appendix C. For instance, the energy of the ground state for 
p > p c3 with the boundary solutions kL, k% and vi, is transformed into that 
of p < —pc3 with k x , k 2 and v v However, it seems quite non-trivial how the 
two sets of the charge rapidities for the two cases of p > and p < could 
be related to each other. (See also Appendix C.) 
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Let us show that momentum Ul which is real- valued when p < p c \ be- 
comes complex- valued when p > p c \. First, we note that when k is real and 
1 7T — k\ <C 1, we have 

1 , / 1 + pe tk \ TT . . 2 psm(n-k) . 

Here H(x) denotes the Heaviside step-function: H(x) = for x < and 
H(x) = 1 for x > 0. Suppose that momentum fc^ be real even when p > 1. 
Since is close to 7r, we have Z L (kL) = /cl/7t + z c B (kL)/L. It follows from 
(H) that the value of z B (ki) for p > 1 is by 1 smaller than that of the 
case when p < 1: z c B {k L ) = k L /n — 1 + 0(1/L) for p > 1. Thus, we have 
/ L /L = k L /ix + (A; l /tt - 1)/L + 0(1/L 2 ), which leads to k L = vr + 0(1/L 2 ) 
for II = L. However, when k = 7r the wave function should vanish under the 
open-boundary condition. Thus, we arrive at an inconsistency. Therefore, 
the momentum k^ should be complex- valued when p > 1. 

We can show that V\ becomes imaginary when p > p c2 . Let us take the 
following branch of the logarithmic function: — ilog (e(x)) = it — 2tan _1 (x), 
where e(x) denotes e(x) = (x + i)/(x — i). Then, we can show 

Y ( v + riyu) = I 6 ^(v) + n -,(v), for 1 < n, 
nl + 7 j \6, +n (v) + l-d^ n (v), for 7 >n. 

Applying the formula (^) with 7 = a to the function z s B (v), we can show 
that if we assume the smallest rapidity v\ to be real, then it would be 0(1/ L 2 ) 
for p > p C 2, and also that therefore it should be imaginary when p > p C 2. In 
the same way with the rapidity v\, using the formula (|T0|) we can also show 
that momentum ki becomes imaginary when p > p c3 . 

We can evaluate the largest and smallest integers of all the possible quan- 
tum numbers for real momenta in the following way. The function Z1(k) is 
monotonically increasing with respect to k, since the density of real momenta 
should be non-negative. We note that under the open boundary condition, 
the Bethe-ansatz wavefunction should vanish if there exists a momentum of 
k = or n. Thus, the equations for I min and I max are given by 

Z C L (0) = (I mm - l)/L, Z c l (tt) = (I max + I)/ 'L. (11) 

We determine I min and I max by solving eqs. ([II]). For instance, for the case 
when < p < 1, it is easy to see Z L (0) = and Z L (n) = (L + l)/L, so that 
we obtain I min = 1 and I max = L. 
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For real rapidities, we can obtain J min and J max by applying the argument 
in Ref [|Tj|]. It is easy to show that they satisfy the following equations. 

Z' L {0) = (Jmin - I)/ L, Z s L {oo) = (J max + 1)/L. (12) 

Solving eqs.([L2|) we determine J m i n and J max - For example, let us consider 
the case p c2 < p < p c3 . From eqs. and (H) we can show Z£(oo) = 
1 + (1 — Jmax) I ' L- Thus, we obtain J max = L/2. We can discuss the maximal 
and minimal quantum numbers also for some excited states with boundary 
solutions, similarly. Some details are given in Appendix D. 

The new hole appears in the half-filled band, when p > p c3 . Therefore, 
there is a gapless mode of particle-hole excitations for the half-filled ground 
state under the open boundary condition. Let us give some explanation in 
three paragraphs in the following 

First, we consider the appearance of the new hole. This is a consequence 
of the formation of the boundary k — A string. In fact, the number of possible 
real momenta in the band is given by L — 1, since I max = L and J m j„ = 2 
when p > p c3 . On the other hand, there are only L — 2 real momenta in the 
wavefunction since we have two complex momenta k\ and ki- Thus, there 
should be one hole in the band. By shifting the quantum number of the hole, 
we can make a series of charge excitations with the hole; for the ground state 
the quantum number of the hole is given by L, while for the charge excited 
state it is given by an integer less than L. 

Second, we consider the gap energy for the charge excited state, where 
the quantum number of the hole is close to L (for example, L — l). Then, we 
see that it becomes infinitesimally small when we take the thermodynamic 
limit L — > oo. In this sense, we may call the mode gapless. Here we note 
that the excitation energy should be continuous with respect to the charge 
rapidity kh of the new hole. 

Furthermore, we can explicitly calculate the charge excitation energy, 
applying the method |T{| of the finite-size correction. We recall that kh 
denotes the charge rapidity of the new hole. We denote by Ef^ikh) the 
energy of the charge excited state with the new hole. Then it is given by 

Ef(k h ) = E{- 2e c (k h ) + 2e c (7r), (13) 

where E 9 L denotes the ground-state energy for p > p c3 and e c {k) is the dressed 
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energy for the half-filled band given by 

A r°° e- ulJJ Ji(^)cos(^sinA;) 

e (k) = cosfc- / — — — -duj. (14) 

2 J-oo cjcoshwu; 

The expression of the chemical potential A at the half-filling will be given 
later in ( |2"4] ) . Q From the expression of the excited energy (|I^) we see that 
the gap energy of the mode is of the order of 1/L 2 . Thus, we see that the 
gap energy vanishes under the thermodynamic limit: L — > oo. 

The three boundary complex solutions for p > p c3 can be considered as a 
variant of k — A string that was originally defined for the periodic Hubbard 
model. In fact, we can derive the expressions (§) of the boundary k — A string 
from the viewpoint of classification of k — A strings of length n — 1. Details 
will be discussed in later papers. 

Let us explicitly study for a finite-size system the behaviors of momenta 
and rapidities with respect to the boundary chemical potential. In Fig. 1, 
the flows of momenta and rapidities are plotted versus the parameter p for 
the 8-sited Hubbard Hamiltonian under the open boundary condition, where 
the roots are obtained numerically by solving the Bethe ansatz equations 
with L = N = 8 and M = 4. As far as the finite-size systems we have 
investigated are concerned, the numerical solutions are consistent with the 
following consequences of the analytic approach: the complex solutions are 
formed one-by-one at the critical points of the parameter p\ there is a charge 
hole when p > p c %. This is nontrivial. The analytic method should be valid 
only when the system size is very large. However, these important properties 
are already observed in such a small system as the case of L — 8. 

Fig.l 

Let us explain how we apply to our system the method |19| of the finite- 
size correction. We consider the Hamiltonian 7i = 7i — AN — h(N — 2M)/2, 
where A and h are the chemical potential and the uniform magnetic field, 
respectively. In order to define densities of the roots of the Bethe ansatz 
equations, we extend Zl(k) and Z s L {y) into odd functions defined both on 
positive and negative values of their variables. For an illustration, we consider 
the density of real-valued rapidities. When < p < p c2 , we have Zf,(0) = 

2 The expression ([l3]) can be derived from the formula (|l^) by replacing the hole 
momentum k g h of the ground state by that of the excited state. 
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and the function Z s L (v) itself can be simply extended into an odd function of 
v by Z s L (—v) = —Z s L [v) for v > 0. We define rapidity with negative suffix by 
v - m = —v m for m = 1, . . . , L/2. Then, the density of the real rapidities is 
given by the derivative p s L (v) = dZ s L {v)/dv for — oo < v < oo. When p c2 < p, 
however, the function does not vanish at the origin: Z s L {ti) = 1/L. In this 
case, we introduce some shifts of the function and the variable Z s L (v) = 
Z s L {v) — 1/L and v m = v m+ i, respectively. We also introduce rapidity of 
negative suffix by V- m = v m , for m > 0. Then, the Bethe ansatz equations 
are given by 

Z S L {v m ) = J m /L, for m = —J^axi ■ ■ ■ > ^maxi (15) 

where J m = m and J^ ax = J max ~ 1- Then, we can safely define the density 
of rapidities by the derivative p s L (v) = dZ s L {y) / dv . 

Taking the derivatives of the Bethe ansatz equations together with some 
continuous limits, we can systematically derive a set of equations for the 
densities of the system with L sites. For the half-filled band under zero 
magnetic field, the set of equations for the densities up to 0{ 1/L) is given in 
the following 

I 1 POO 

pF L (k) = - + -T c{0) (k) + cosk a 1 (smk-v)p s L (v)dv + 0(l/L 2 ), 

7T L J-oo 

p s L (v) = ^t s{0 \v) + ai (v -sink)p c L (k)dk 

a 2 (v-v')p s L {v')dv' + 0{1/L 2 ). 

-oo 

(16) 

Here a n (x) is defined by 2na n {x) = (2nu)/(x 2 + (nu) 2 ). The boundary terms 
T c (°\k) and t s ^(v) in eqs. ( |i6|) are given by the derivatives of Po(/c)/(27r) 
and Qo(v)/(2ir), respectively, where they are related to z c B {k) and z s B {v) by 

P (k)/2n = z c B (k) - ^(sinfc), Q (v)/27r = z s B (v) - 9 1 (y) + 6 2 (v). (17) 

We now evaluate the ground-state energy E 9 L of the Hamiltonian TC at 
half-filling under zero magnetic field. From eqs. (^) we have the following 

E 9 L = - 2 cos kj - 2 cos k i - AN 

= L eoo + l + A/2 + (e,r^)- £ 2e c (k 9 h ) 

hole 



10 



- (A + 2 cos kj) + 0(1/ L), (18) 

j'eA c 

where T^°'(k,v) = (r c( -°\k), t s (°\v)) denotes the surface density, the symbol 
e = (e c (fc), e s {y )) denotes the dressed energy |L9], |ll| We recall that fc^'s 
denote the momenta of possible holes at the ground state. (For the ground 
state of p > p C 3, we have only one hole. ) The inner product (e, t^) is 
defined by the following [|l^, O 



poo 

e c (k)T c(0 \k)dk+ / e s (v)T s(0) (v)dv. (19) 



Let us define the surface energy e sur of the system by the 0(1) part of the 
ground-state energy. Then it is given in the following 



1. For < p < 1, 

n=0 



~ 2 n /• 0O 2e-^J 1 MJ 2n (w) , 

econ + p-Y^P V —doo. 20 

Jo uj cosh mu; 



2. For 1 < p < p c3 , 

2e-^cosh(a;sinh/c)Ji(u;) ~ 1 2e- UUJ J 1 (cu)J 2 n(uj) 



r°° 2e "~ cosh u; smh « Ji w , ~ 1 /■ 

e co „+p- / ^-r ^— — / 

jo u; cosh MtU rZi P Jo 

3. For p c3 < p, 



n=l P 2U J° w cosri WCJ 

(21) 



du. 



1 « 1 «o 2e-" w Ji(o;) J 2n (o;) , 

e con + 4 M -A-- + ]T— / u ; 2nl (22) 

P n^i V jo cosh wo; 

Here the symbol e con denotes the surface energy for p = 0, which is explicitly 
given by 

, roc p -2uw j ( \ 

e con = (1 - A/2) + 2 V / TT^ - 2u - / -i^da;. (23) 

Jo a; cosh uu 

The chemical potential A at half-filling is given by 

^2-2^1'^. (24) 
Jo cosh ULO 
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Let us discuss the ground-state energy for the strong-coupling case. When 
p > p C 3, it becomes close to the energy of the first charge-excited state for 
p = 0. We compare the surface energy for p = given by eq. (^) with that 
of p > p c3 given by eq. (0). Then, the main part of the difference between 
them is given by 4u = U, which is almost equivalent to the charge-gap energy 
Au — 2A at p = 0. We note that when u 1, we have p c3 ^> 1 and u ^> A. 

Under the strong coupling condition, the main part of the surface energy 
is given by the following; p + e con when 1 <C p < p c3 and 4u — 2 A + e con when 
V > Pc3- /,From the calculation of the ground-state energy, we can evaluate 
the average number n% of electrons on the 1st site, since it is defined by ri\ = 
dE L /dp. We find that dE L /dp m 1 for 1 p < p c3 and dE L /dp for p c3 < 
p. This suggests that one hole should be localized at the surface site when p > 
p c3 . The result is consistent with the discussion over the complex boundary 
solutions that the half-filled ground state has gapless charge-excitations when 
p > p c3 since one hole appears in the band. 

Let us discuss the spectrum of a finite-size system numerically. The low- 
excited spectrum of the 6-sited open Hubbard Hamiltonian with U = 20t is 
obtained by the exact numerical diagonalization of the Householder method. 
The spectral flows with respect to the parameter p are depicted in Fig. 2. 

Fig. 2 

From Fig. 2 we see that the energy levels of charge excitations become close 
to the ground-state energy at p = p c3 . 

From Figs. 1 and 2, we have the following observations, respectively, 
(i) When p > p c3 the first charge-excited state can be obtained by shifting 
the position of the hole in the band of real momenta k/s; such shifting is 
equivalent to taking a different quantum number for the hole, (ii) The energy 
level of the first charge-excited state for p > p c3 is identified with that of the 
lowest state above the charge gap for p = 0; we can trace the spectral flow of 
the excited state from p = p c3 down to p = in Fig. 2. ^From the analytic 
approach, the observation (ii) should hold due to the adiabatic hypothesis on 
the quantum numbers. From (i) and (ii), we can say that the characteristic 
properties of the energy spectrum discussed by the analytic method are also 
in common with that of the finite-size system. Thus, the spectrum of the 
6-sited system may illustrate that of thermodynamically large systems. 

In this paper we have discussed the boundary solutions for the half-filled 
band when p > 0. We have shown that when p > p c3 one mode of charge 
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excitations has the gap energy of the order of l/L 2 ; we call it massless since 
the gap vanishes in the thermodynamic limit. By the method of the finite- 
size correction, we have calculated the ground-state energy and the excited 
energy of the massless mode. We note that it is not difficult to derive explicit 
formulas for the energies of other excitations. In fact, all the spectrum shown 
in Fig. 2 can be explained analytically. Details will be given in the next 
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Figure Captions 



fig.l 

(a) Flow of momenta (or charge rapidities) as a function of boundary po- 
tential {p) for the one-dimensional open Hubbard model with 8 sites at the 
half filling (N = L) and U = 201 Solid lines represent real momenta. At 
P = Pel ~ t, the largest charge rapidity approaches n and becomes complex, 
7r + i/t, whose complex part is given by the dashed line, for p > p cl . Beyond 
P — Pc3 ~ U, the smallest momentum becomes complex, —in', where k' is 
represented by the dot-dashed line, (b) Flow of rapidities (or spin rapidities) 
for the same system. Solid lines represent real rapidities. At p — p c2 ~ U/2, 
the smallest rapidity becomes complex, i%, where x is given by the dashed 
line. 

Fig.2 

Spectral flow for the 6-site open Hubbard chain at the half-filling with U = 
20t as a function of boundary potential p. Dots denote all of the eigenvalues 
for this system obtained by the direct diagonalization. The lower solid line 
represents the ground state energy given by the Bethe ansatz. The upper 
solid line corresponds to the first charge-excited state, which can be traced 
back from p = p c3 to p = 0; at p = 0, it is the lowest level beyond the charge 
gap. Enlarged flow around the gap-closing transition point is depicted in the 
inset. Below a critical point (p^), a charge gap exists above the continuum 
of low-energy spin excitations. 
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1 Appendix A: 



Let us briefly outline the derivation of the Bethe-ansatz equations through 
the algebraic Bethe-ansatz method for the open-boundary XXZ model given 
by E.K. Sklyanin. Some details can be found in Ref. pTl[] . (See also |13| .) 



We write the eigenstates for N electrons with M down-spins as 

^nm = x N )cl iai ■ ■ ■ c\ NON \vac). (A.l) 

Here, Xj and <7j are the position and spin variables of the electrons, respec- 
tively. In the region xqi < • • • < xqn, we assume that the Bethe-ansatz 
wavefunction / takes the form 

N 

/<T1,...,(TJV (x 1} - ■ ■ ,x N ) =J2 e P A a Q1 ,-,a QN {kpi, ■ • • , k PN ) exp{i k Pj X Qj}- 

P j=l 

(A.2) 

Here the Q is an element of Sn, the permutation group of iV particles, and P 
runs over all the permutations and the ways of negations of k's; there are TV! x 
2 N possibilities for P, while N\ for Q. We employ the notation: k_j = —kj. 
The symbol ep denotes the sign of P; if the permutation is even, P makes 
ep = —1 when odd number of fc's are negative and ep = 1 when even number 
of fc's are negative. Let us introduce the vector A(kj 1 , . . . , kj N ) such that its 
element for entry (Ql, . . . , QN) is given by A aQlt ... <aQN (kj 1 , • • • , kj N ). Here we 
note that the suffix ji, . . . ,jn can be written as PI, ... , PN, respectively, by 
some P. Then, we can show that the consistency condition for the amplitudes 
A(kj 1 ,. . . , kj N ) is given by the following 

T(smkpi)A(k P1 , . . . ,k PN ) = A(k P1 , . . . , k PN ). (A. 3) 

Here T(u) is the inhomogeneous transfer matrix of the open-boundary XXX 
model with inhomogeneous parameters sinfcpi, . . . , sin kpN- Let us 
denote the eigenvalue of the transfer matrix T(u) by A(u). Then, from the 
condition A (sin kpi) = 1 and the Bethe ansatz equations for the XXX model, 
the Bethe ansatz equations of the ID open-boundary Hubbard model for the 
charge and spin parts are obtained, respectively. 
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2 Appendix B: Stability of the boundary so- 
lutions 

Let us discuss explicitly the stability of the boundary solutions appearing 
in the ground state and some excited states, both for the cases p > and 
p < 0. We recall some notation in the following. We have defined the symbols 
k and ot by n = \p\ and a = sinh k/u, respectively. The symbol p C j is given 



by p C j — (j ' — 1) M + y 1 + (j — l) 2 u 2 for some integer j; we note that p = p C j 
corresponds to a = j — 1. 

Let us introduce a useful formula in the following 

i I (ry _|_ 77^2^2 _|_ ^2 \ 

n (i7W-v) + n (i7it + v) = — In — - , for v > 0, 

27r y (7 — n) z u z + v l ) 

(B.l) 

where 7 > 0. We can show eq. ( |B.1|) by a similar method for the formula 
(|T0D; we take the branch of the logarithmic function, and use the relations 
arg{i{pi — n)u — v ) = arg{v — 2(7 — n)u) — n and arg{i{pf + n)u — v ) = 
arg(v — 2(7 + n)) + ir. 

For the case of p > 0, we may consider the three complex roots k\, ki 
and Vi in the following regions. 

ki = 7T + ik — i5l, for a > 0, 
t>i = i(o; — l)u + zr/i, for a > 1, 

fci = z log f(a - 2)u + y/(a - 2) 2 v? + 1 j + ft, for a > 2. (B.2) 

We call the boundary solutions stable when ft, 8\ and 771 are very small. For 
some convenience, we use symbols e\ and defined in the following 

sin &i = i(a — l)u + ie\, sin k L = i(a — l)u + ie^, (B.3) 

which are related to ft and ft by 



6l = — cosh k x ft, ei = y (a — 2) 2 v? + 1 x ft. 

Let us give explicitly some evaluations of 6l, r]i and e\ for the case of 
p > 0. We assume that p is not close to any of the critical points p c /s. Then 
for the ground state and some excited states we can show the following. 
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1. When there is only one boundary solution ki, we have 

e L = (p- 2L ) . (B.4) 



2. When there are two boundary solutions k^ and v\ and when u > 1, we 
have 



a -~ 2N " 



2- a 



a ^ 2N 



3. When there are three boundary solutions ki,ki and v\ and when u > 1, 
we have 





-ml 


= o( 


kll 1 




-ml 


- o| 


N 




£l| 


- o| 


N 



2L 



a 



1 



2L / a \ -2L 



2 A 



p- 2L ■ (B.6) 



Here denotes the following 



zi = exp(ifci) = -(a - 2)u + y/(a - 2) 2 u 2 + 1. 
It is easy to see that \z%\ < 1 for a > 2. 

We note that the evaluations ( p.4| ), (p.5| ) and (p.6| ) can be applied for the 
half-filled ground-state solutions in the regions of p c i < p < p& (0 < a < 1), 
Pc2 < P < Pci (1 < a < 2) , and p& < p (2 < a), respectively. 

In the derivation of (|B.5|) and ( B.6|) , we have assumed that u > 1. In 



fact, applying the formula ( |B.1| ) with 7 = a — 1, we can make the following 
approximation when u > 1 



/ , / a 2 u 2 + sin 2 kj 
, ,_±i ; 2tt"" V(a-2)% 2 + sin 2 fc j< 



53 53 ^(*(«- l)«-rsin%) = 53 "T ln I WvT77i ^2 



2^ ln • (R7) 



18 



We recall here that N is the number of electrons. The approximation ( |B.7|) 
can be not effective when u is very small. However, it seems that it is 
nontrivial to evaluate and rjx for the weak-coupling case; more precise 
estimates on sin 2 k/s should be necessary when u is very small. 

For the case when ( |B.5| ) is valid, the quantity el should be very small if 
the following inequality holds 



V 



_1 OL 



2 - a 



< 1. (B.8) 



For the region: p c2 < p < p C 3 (1 < a < 2) , however, the inequality ( |B.8|) does 
not hold for all values of a satisfying 1 < a < 2. Let us consider the case of 
u ^> 1; when u is very large, we can approximate 1/p by l/(2cra) using the 
relation: p = au + a/1 + a 2 u 2 . Then, we can show that the inequality ( |B.8| ) 
holds under the following condition 

a<2 . (B.9) 

2u 

Thus, at least for the case of u ^> 1, we have shown that the boundary 
solutions k L and V\ are stable when 1 < a < 2 — l/(2u), where 2 — l/(2u) 
is very close to the critical point a = 2. For the region: p > p c3 (a > 2), we 
can show, under the condition: u ^> 1, that the boundary solutions and 
v\ are stable if a > 2 + 1/ (2w). 

Similarly, we can discuss the case when the evaluation (|B.6|) is valid, 
where there are the three boundary solutions, ki, k\ and v±. For the strong- 
coupling case, we can explicitly show that the boundary solutions are stable 
if a > 2 + l/(2w); we have the following 

p -1 ^! -1 < 1, when a>2 + l/(2u). (B.10) 

a — 2 

Let us now consider the boundary solutions for the case when p < 0. 

k 1 = in — i5 1 , for < a < 1, 

v 1 = i(a — l)u + ir)i, for 1 < a < 2, 



fc 2 = ihg((a-2)u + ^/(a-2) 2 u 2 + 1) +i5 2 , for 2 < a. (B.ll) 
Here 5 2 and should be very small. We recall that for the case of p < 0, 



the boundary solutions k l} k 2 and v 1 have been discussed in Ref. |TjJ. For 
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some convenience, we use symbols t x and e 2 denned in the following 

sin k x = i(a — l)u + ze 1; sin k L = i(a — l)u + ie 2 , (B.12) 
which are related to 5[ and S 2 by 



e x = — cosh k x 5 l} e 2 = y (a — 2) 2 u 2 + 1 x 5 2 . 

Applying the formula (|B.1|) , we can evaluate e 1; r} x and e 2 as follows. 
1. When there is only one boundary solution k' x , we have 

e[ = O (p- 2L ) . (B.13) 



2. When there are two boundary solutions k[ and v[ and when u > 1, we 
have 



a -- 2N ' 



2- a 



a ^~ 2N 



O ( ( ^f - ) p-" ] . (B.14) 



3. When there are three boundary solutions k[, k' 2 and v\ and when u > 1, 
we have 



o(i4l 2L ), 



Iq-^il = O k: 



' |2L 



ft 



a - 2 



-2iV N 



"ii - (B15) 



Here z 2 is given by 

z 2 = exp(ik 2 ) = —(a — 2)u + \J (a — 2) 2 u 2 + 1. 
We note that \z 2 \ < 1 when a > 2. 
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We note that the evaluations ( |B.13| ), ( |B.14| ) and ( |B.15| ) can be applied for 
the ground-state solutions in the regions of —p C 2 < p < —p c \ (0 < a < 1), 
—Pc3 < p < —Pc2 (1 < a < 2), and p < — p c3 (2 < ft). We also note that in 
the derivations of ( |B.14|) and ( |B.15|) we have made the same approximation 
with ( |B7D . 

For the case when the evaluations (p,13| ), ( p. 14]) and (P-15| ) are valid, 



the stability conditions for the boundary roots are satisfied for any p and u. 
For the case we have the following relations when p < — 1 and u > 



br 1 < i, [p- 1 - 



< i, \p x 4 



ft -2 



< 1, etc. 



(B.16) 



However, we should remark that it is not certain whether (|B. 14|) and (p,15| ) 
are valid also for the weak-coupling case: u<l. 

We now discuss how the boundary solutions of the open-boundary Hub- 
bard model can be related to those of the interacting spin- 1/2 fermion sys- 
tems. We consider the case when the band width t is very large and the 
electron density N/ L is very small. 

In order to show explicitly the effect of the large band-width, we replace 
u and p in eqs. (H), © and (|5[) by u/t and p/t, respectively. We recall that 
so far the energy scale has been normalized such that t = 1. Under the limit 
of t — > oo, the critical points p' c s become the following 



Pcl/t 
Pa/t 



1+u/t, 
1 + 2u/t. 



(B.17) 



The values obtained in the limit are equivalent to the critical points of the 
boundary parameter [17] for the interacting spin- 1/2 fermions. 

When the electron density N/L is very small, the Fermi wavenumber Jcf 
should be very small. Therefore, we can make linear approximations for the 
charge rapidities such as exp(ik) ~ 1 + ik and sin/c w k. Then, we can 
show that the boundary term in the Bethe ansatz equations of the open- 
boundary Hubbard model corresponds to that of the interacting spin-1/2 
fermion system under the linear approximations 



2k 1 / 1 + exp(ik)p/t \ 
2tt 2ni 11 \1 + exp(~ik)p/t ) 



2tt 



tan 



-(i+p/0 



(B.18) 
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Hereafter we may renormalize the boundary chemical potential p so that we 
can replace p/t by p. Thus, we have explicitly shown that when f > 1 and 
N/L <C 1, the Bethe ansatz equations of the open-boundary ID Hubbard 



model are reduced into those of the interacting spin-1/2 fermion system |L7 
with the open-boundary condition. 

Under the large band- width and small electron-density limit, the bound- 
ary solutions of the open-boundary Hubbard model for the case of p < re- 
main intact. We can make the same approximation with (p. 7f ), since sin 2 k/s 
are very small. Here we note that the case of large band-width corresponds to 
the weak-coupling case where the approximation ( |B.7| ) can be non-effective. 
When the density is very low, however, it is valid for some cases. For ex- 
ample, we may consider the case where N is fixed and L is proportional to 
t under the limit t — > oo so that each momentum kj is proportional to 1/t. 
Then, we can apply the approximation ( |B.7| ) for this case. In this way, the 
boundary solutions of the open Hubbard model for the case p < are related 
to those of the interacting spin-1/2 fermion system discussed in Ref. |17 |. 



For the case when p > 0, however, the boundary solutions of the open 
Hubbard model are not related to any solution of the interacting spin-1/2 
fermion system. They exist only when the band is half-filled. The physical 
condition is completely different from the low density case. 



3 Appendix C: Particle- hole transformation 
for the open-boundary Hubbard chain 

Let us denote by dj >a and o?] )fT , the annihilation and creation operators for 
the hole with spin a on the jth site, respectively. We define a particle- 
hole transformation by the following. We replace the creation (annihilation) 
operator of electron with spin a on the jth site by the annihilation (creation) 
operator of hole with spin a on the jth site for cr =|, J, and for j = 1, . . . , L, 
and then multiplying the gauge factor (— 1) J to the hole operators on the jth 
site for over all the sites: 

4 - (-i)^>, CijCT - (-i)xV (c.i) 

The ground-state energy for p > is related to that of p < by the 
particle-hole transformation; the sign of the boundary chemical potential 
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is changed under the transformation. Let us denote by E(N±, Nf,U,p) the 
ground-state energy for down-spin electrons, iVj up-spin electrons with 
the Hubbard coupling U and the boundary chemical potential p. Then, 
applying the particle-hole transformation , we have the following 

E(L- M,L- M';U,p) = E(M, M' ';[/, -p) + (L - N)U + 2p. (C.2) 

For the half-filled band, the ground-state energies for p > and p < are 
explicitly related. Recall we assume L is even. Then, we have the following 

E(L/2, L/2; U, p) = E(L/2, L/2; U, -p) + 2p. (C.3) 

On the other hand, it seems quite difficult to find out an explicit relation 
between the sets of the charge (spin) rapidities for the cases p > and 
p < 0. It seems as if there might be such a simple connection that for any 
momentum k in the ground state of p > the value is ± k corresponds to one 
of the ground-state solutions for p < 0. However, it is not the case. There is 
no such relation between the boundary solutions k ± , k 2 and v 1 for p < — p c3 
and the boundary solutions k±, k^ and V\ for p > p c3 . 

From some numerical solutions of the Bethe ansatz equations, it is sug- 
gested that it can be quite nontrivial to find out any explicit relations between 
the sets of the half-filled ground-state solutions for the cases p > and p < 0. 
Some details should be discussed in later papers. 

4 Appendix D: Boundary solutions for some 
excited states 

We discuss the quantum numbers of some excited states with the boundary 
solutions. We assume the adiabatic hypothesis for the quantum numbers. 

Let us consider an excited state which have only real-valued momenta 
and rapidities when p = 0. We denote by Aq ( Aq) the set of the quantum 
numbers for the momenta ( rapidities). Let us denote by A C (A S ) the set of 
the quantum numbers of real-valued momenta in the excited state at a given 
value of p. In general, A C (A S ) depends on p. For the excited state we write 
by A^ m (A| m ) the set of the quantum numbers of complex-valued momenta 
(rapidities). It is useful to introduce the notation for holes; we denote by 
A c hole (A s hole ) the set of the quantum numbers of holes for the real-valued 
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momenta (rapidities). Then, we can define Z^(k), Z[(v), z B (k), and z B (v) 
also for the excited state; in the formulas (§) we replace A^ and A^ by A c and 
A s , respectively. Similarly to the ground state, we can evaluate the I m j n 's 
for the excited state as follows. 

Imin — Z B (0) + 1, I m ax = L + Z B {7T) — 1, 

Jmin = z b(®) 1j 

Jmax = (N-N im )-(M-M im ) + (z s B (oo)-l/2). (D.l) 

Here N im and M im denote the number of complex-valued charge and spin 
rapidities (boundary solutions), respectively. We recall that N and M denote 
the number of electrons and that of down-spins, respectively. 

For an illustration, let us discuss the boundary solutions of an excited 
state for the case when p < 0. Hereafter we assume N < L. We consider the 
excited state of N electrons with M down-spins where the quantum number 
at p = is given by the following. 

A£ = {1,3,4,. ..,N + 1}, A* = {1,2,...,M}, 
A im = AL = 0- (D.2) 

It follows from (p.l|) that when p = the sets of holes are given by 

A^ e = {2,iV + 2,iV + 3,...,L}, A s hole = {M + l,...,N-M}. (D.3) 

Here we note that when L = N + 1, then we have A£ oZe = {2}, and also 
that when N is even and M = N/2, then we have A^ o/e = cf). Applying the 
formulas ( |10D and ( p . 1|) , we can show that there are four critical points given 
by —Pcj for j = 1, . . . , 4. We have the following five cases when p < 0. 

1. For — p c i < p < 0, we have no boundary solution. We have a hole at 
1 = 2. 

A c = {l,3,4,...,iV + l}, A* = {1,2,...,M}, 
A^ e = {2,iV + 2,...,L}, A s hole = {M+l,...,N-M}. 

2. For —p C 2 < p < — Pel, we have k[ and a hole at J = 2. 

A im = {1}, A im = 4>, 
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A c = {3,4,...,iV + l}, A S = {2,3,...,M}, 
Al ole = {2,N + 2,...,L}, A s hole = {M + l,...,N-M}. 

3. For —pa < p < —Pc2-, we have k\ and v[ and a hole at I = 2. 

AL = {1}, AL = {1}, 

A c = {3,4,...,iV+l}, A* = {2,3,...,M}, 
A c hole = {2,N + 2,...,L}, A s hole = {M + l,...,N-M}. 

4. For —pd < p < — pa, we have k[ and v[ but no hole at / = 2. A new 
hole appears at / = L + 1. 

AL = {1}, AL = {1}, 

A c = {3,4,...,iV + l}, A S = {2,3,...,M}, 
A c hole = {N + 2,...,L + l}, A s hole = {M + l,...,N-M}. 

5. For p < —pci-, we have k\ and v[. A new hole appears at J = 1. 

AL = {1}, A| m = {1}, 

A c = {3,4,...,iV+l}, A S = {2,3,...,M}, 
A c hole = {N + 2,...,L + l}, A s hole = {l,M + l,...,N-M}. 



We recall p c4 = 3m + \Jl + (3k) 2 . 

Let us consider the ground state and the excited state discussed in the 
last paragraph. For the two regions — p ci < p < —p C 3 and p < —p c ^ the 
ground-state solutions have the same structure, while the excited state have 
the different structures; the excited state has the two boundary solutions 
both for the two regions, however, it has the different numbers of holes in 
the spin rapidities for the two regions. Thus, it is suggested that there can 
be more subtle points in the boundary excitations than had been described 
in Ref. [0 for those of the interacting spin-1/2 fermion system. However, it 
seems that some physical interpretations similar to those in Ref. jr7| should 
be valid also for the boundary excitations of the open-boundary Hubbard 
model. Some precise investigations should be discussed in later publications. 
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